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Abstract 

We investigate the possibilities to calculate vector partition functions by 
means of iterated partial fraction decomposition, as suggested by Beck 
(2004). Particularly, for an important type of families of rational func- 
tions, we describe an algorithm to compute the numerators in their partial 
fraction decomposition "formally," i.e., as a formal expression in the pa- 
rameter. We also analyze the type of generalized rational functions that 
appear during the execution of an algorithm based on iterated partial 
fraction decomposition and explain how to handle these objects. 

1 Introduction 

We denote the set of nonnegative integers by No and the set of nonnegative 
real numbers by R + . Let ci,.. . ,c<j S Z m C R m , contained in an open half- 
space. The typical situation to have in mind is that d > m, so that the family 
(cj) is linearly dependent. To this data, associate the vector partition function 
<p Cl ,.., <Cd : Z m ^ N by 

<Pci,...,c d ( b ) = \i X E N : x l c l "I 1" x dCd = b}\. 

So fci,...,c d {v) is the number of decompositions of b as a nonnegative integral 
linear combination of c\, . . . , c^. In a more geometric language, <p Cl ,...,c d (b) is 
the number of integral points in the parametric polytope 

P Cl ,...,c d ( b ) = {x £ R+ : xxa H h x d c d = b}. 

Thus the theory of vector partition functions generalizes Ehrhart theory [TCJj . 
the latter corresponding to the study of the restriction of (f ci ,...,c d to a ray. 

Vector partition functions occur throughout mathematics as a convenient 
way to encode polyhedral models, such as: 1. In representation theory, as 
Kostant's partition function [T2]. 2. Also in representation theory, as the vector 
partition functions constructed by Billey et al. |5J from Gelfand-Tsetlin pat- 
terns and by myself [7 from Littelmann patterns and Berenstein-Zelevinsky 
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polytopes. 3. For counting integer flows in networks, as discussed by Baldoni et 
al. [2]. 4. In the optimization of computer programs (see |17j for an extensive 
list of applications in this area) . 

The usefulness of expressing a given function by a formula involving vector 
partition functions partially relies on the following theorem, proved by Black- 
ley [B], Dahmen and Micchelli [S], and Sturmfels [15]: Given ci,...,a, there 
is a fan F in R d such that the restriction of (p Cl ....,c d to each maximal cone 
C e F is given by a quasi-polynomial fc, and if Cl ,...,c d vanishes outside F. 
Suppose you are given a formula for some interesting quantity involving a vec- 
tor partition function. If you succeed in determining explicitly its fan F and 
its quasi-polynomials (/c), you obtain a very explicit general formula for the 
original quantity. For this approach to work, it is of course essential to have 
effective means to determine F and (fc)- To date, there are a number of al- 
gorithms known allowing to do this, employing a variety of ideas. Some of the 
approaches I am aware of are the following: 

1. Using Barvinok's algorithm [3]- This has been worked out by Verdoolaege 
et al. PI]. 

2. Using the Szenes-Vergne residue formula [TB]- This was employed by Bal- 
doni et al. [TJ for Kostant's partition function associated with classical 
root systems and by myself [7] for vector partitions functions obtained 
from generalized Gelfand-Tsetlin patterns. 

3. Using Paradan's wall crossing formula which relates the quasi-polynomials 
of adjacent maximal cones of F, see Boysal and Vergne [5]. 

4. Iteratively rewriting the generating function as a sum of rational functions 
with a special pole structure. Milev is working on this approach [14J . Note 
that he uses the term "partial fraction decomposition" in a different sense 
than we do. 

Beck suggested a further algorithm [4], based on iterated partial fraction de- 
composition of the generating function. The aim of this extended abstract is to 
report on work in progress to work this idea out. 



2 Preliminaries 

This section is a summary of Beck's suggestions [4 to compute vector partition 
functions. We denote the set of positive integers by N. Let m, d € N and let 
A G Z( m - n ) be an (m x <i)-matrix with integral coefficients. We consider A as a 
linear map R d — > R m . Suppose that ker(A)nR+ = {0}. Then we can associate 
with A the vector partition function ipA ■ Z m — > No by 

^ A (b) = \{xeN d :Ax = b}\. 

Let Ci , . . . , Cd be the column vectors of A. Then ifA as defined above coincides 



with <p ci ,...ca as defined in § 1 We use the standard multiexponent notation, 
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z b = z bl ■ ■ ■ zfyf. The generating function of tp A , 

f A (z) = <PA(b)z b , 

converges on {z G C m : \z Ck \ < 1 for k = 1, . . . , d} to the rational function 

By definition of the generating function, <£>a(&) is the z^-cocfficicnt of f A (z). 
Equivalently, it is the constant coefficient of f A (z)z~ b , 

tp A {b) = const f A (z)z~ b . 

Consider f A (z)z~ b as a rational function in one variable z m over the field 
K m _i = C(zi, . . . , z m -i). Assume there is a well-defined notion of the constant 
coefficient of a rational function. Denote the constant coefficient of f A (z)z~ b G 
K m -i(z m ) by const Zre f A (z)z~ b G K m _\. Then by iterating this procedure we 
get the expression 

</?a(&) = const 2l • • • const Zm f A (z)z~ b G ifo = C 

for the vector partition function ip A . Hence, the problem of computing Lp A is 
reduced to the problem of finding the constant coefficient of a given rational 
function in one variable. This can be done using partial fraction decomposition, 
as described in the following. 

Let if be a field. Consider the field K(w) of rational functions in one variable 
w over K. (In the setting of the previous paragraph, K = C(z\, . . . , z m _i) and 
w = z m .) Let f(w) = G K for polynomials r(w),s(w) G K[w]. Suppose 

deg(r(w)) < deg(s(w)). We want to compute const™ f(w)w~ b as a function 
of b G Z. We can suppose s(0) ^ by a shift in b. Fix a decomposition 
s(w) = 8i(w) • • • Sd(w) into pairwisc coprimc factors. Then Sk(w) and w are 
also coprime for k = l,...,d. Hence, for all b > 0, we get a partial fraction 
decomposition 

f{w)w - b = M^l + ... + ^ + ?<M (2) 
Si(W) s d (w) w" 

for some A k (b;w),B(b;w) G K[w\ such that deg(A k (b; w)) < deg(sk(w)) and 
deg(B(6; to)) < 6. Then 



, x -6 . (&;«>) . . . 4f(&;«>) . , B(b;w) 

const™ j(w)w = const™ — !-••• + const™ — V const™ s — 

s\[w) Sd(w) w" 



Aijb-,0) A d (b;0) 
*i(0) ■" s d (0) 
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Note that, in particular, the constant coefficient of f(w)w~ h does not depend 
on B(b; w) at all, and neither on the higher coefficients of Ak(b; w). Hence we 
do not need to know the complete partial fraction decomposition of f(w)w~ b 
to proceed with our computations. 

By the above ideas it is indeed possible to compute quasi-polynomial ex- 
pressions for vector partition functions. Beck demonstrates this by explicitly 
calculating the case A = ( x x o l ) • However, this example is special in some re- 
spect: 1. The last row of A only contains the numbers and 1. Hence, all roots 
of the polynomials 1 — z Ck S C(zi)[z 2 ] will be elements of C(zi). In general, they 
would be elements of an algebraic extension. 2. All components of A are nonneg- 
ative. Given a matrix A such that if a is defined, one can always suppose that 
all components of A are nonnegative by a unimodular change of coordinates, so 
this assumption is without loss of generality. Still it is possible to define ipA for 
matrices A with negative components, so it would be nice to understand how to 
deal with those directly. 3. The columns of A are pairwise linearly independent. 
This translates to the fact that the generating function has no multiple poles 
when considered as a function in zi . Additional complications will arise if there 
are multiple poles. 

There is further study needed to make the proposed procedure applicable 
mechanically In particular, the following questions must be addressed: 

1. Which decomposition of s(w) should one use? There is a natural one 
coming from 0, and the usual fine decomposition into prime powers. 

2. What is the constant term of a rational function? In general this de- 
pends on the region of convergence considered. If handled without the 
appropriate care, this leads to paradoxical formulas such as 

1 = const (l + z + z 2 + ■ ■ • ) = const 

v ' 1 — z 

— const — z~ x j = const(— — z~ 2 — •■■)= 0. 

1 — z 1 

3. Which algorithm should be used to compute the partial fraction decompo- 
sition? Note that it is crucial to be able not only to compute a particular 
partial fraction decomposition of f(w)w~ b for a a particular b e No, but 
as a function of the variable b. 

4. How do we deal with multiple poles of /^(w)? These occur if and only if A 
has a pair of linearly dependent columns. Hence they become increasingly 
unlikely for higher m, but cannot be excluded. 

5. If we answer the first question by using the fine decomposition of s(w) , lots 
of summands and roots of unity will appear. Do we have to compute all of 
them individually, or can we use the Galois invariance, i.e., the invariance 
of (fiA(b) under the G(Q ab |Q)-operation, to facilitate this? 

We will partially address these questions in the following. The remaining issues 
are delayed to a forthcoming full version of this extended abstract, where the 
word "Towards" will be hopefully dropped from the title. 
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3 A lemma about partial fraction decomposition 



By partial fraction decomposition I mean the following: 

Definition 1. Let O be a domain. Let ^ be an element of the quotient field of 
O. Fix a decomposition s — s\ ■ ■ ■ Sd of s such that (sj, Sk) — O for all j =/= k. 
Then a partial fraction decomposition of - with respect to the decomposition 
si, . . . , sj, is an expression 

- = 1 1 \-p. (3) 

s Si s d 

with n,. .. ,r d ,p £ O. 

Note that in this generality the partial fraction decomposition is not unique. 
Of course, in the important case when O is Euclidean, one can standardize fur- 
ther to obtain uniqueness. Partial fraction decompositions always exist. This 
follows from the following lemma, which gives more importantly a way to de- 
termine the numerators r k explicitly. Note that, once the numerators are deter- 
mined, the integral part p can be obtained by subtracting the fractional parts 
from -. 



Lemma 1. Let r, s, sx, ■ ■ ■ , Sd & O as in Definition 1 Let Ck = — for k £ 
{l,...,n}. Then Ck is invertible in 0/(skf. Let dk G O be an inverse of Ck 
modulo (sfe) and rk = rdk- Then the give rise to a partial fraction decompo- 
sition of i.e., (|3| holds for some p G O. 

Proof. Let k £ As the elements Sj and Sk were supposed to be 

coprimc for j =/= k, these sj are invertible modulo (s/c)- Hence their product Ck is 
also invertible modulo (s^), and we can define and as in the announcement. 
By definition, d^Ck = 1 mod (sfc). Hence 

r k c k = rd k c k = r mod (s k ). (4) 

For j ^ k we have s k \ Cj, so 

rjCj = mod (sfc). (5) 

Combining Q and we get 

T\C\ H V r d c d = r mod (s fc ). 

This congruence holds for all k, so by the Chinese remainder theorem it also 
holds modulo (s). In other words, 

r = r 1 c 1 -\ \-r d c d +ps (6) 

for some p € O. Dividing ^ by s we obtain the announced partial fraction 
decomposition. □ 



5 



Note that Lemma 1 allows for the computation of individual numerators 
Ti without necessity to compute the complete partial fraction decomposition. 
It is also implicitly underlying Kung and Tong's fast algorithm [13 for partial 
fraction decomposition and Xin's algorithm [TS] for MacMahon's partition anal- 
ysis. In my opinion, |LemmaT| is largely underestimated and should, indeed, be 
taught to students as a standard method for partial fraction decomposition. 

We continue by studying the important special case that O — K[w] is a 
polynomial ring and we have a partial fraction denominator of the form 1 — aw n . 
This is the kind of situation typically arising during the computation of a vector 
partition function. 

Definition 2. Let n £ N. Let K be a field of characteristic not dividing n, 
containing the n-th roots of unity. Let a £ (K x ) n be an n-th power. Let 
f(w) € K[w] be a polynomial such that f(a~ 1 ) ^ if a n — a. With these we 
associate the generalized Dedekind sum Sf^ a (n) £ K given by 

Sf.a(n) = - "52 ~F7~ TT • 

a£K J y ' 
a n —a 

Note that Sf ta (n) is the type of generalized Dedekind sums studied by Gessel 



[11 . In this situation, Lemma 1 specializes as follows 



Lemma 2. Let n G N. Let K be a field of characteristic not dividing n, 
containing the n-th roots of unity. Let a G (K x ) n be an n-th power. Let 
f(w) £ K[w] be coprime to 1 — aw". Consider the partial fraction decomposition 

1 A(w) B(w) 



(1 - aw n )f{w) 1 - aw n f(w) 

with deg(A(w)) < n and deg(B(w)) < deg(f(w)). Then the constant term of 
A{w) is A(0) = S f , a (n). 

Proof. We calculate the partial fraction decomposition with respect to the finer 
decomposition 

(1 - aw n )f(w) = f{w) J] (1-aw) 

a n =a 

Let A a denote the numerator of 1 — aw in this finer partial fraction decomposi- 
tion. As deg(l — aw) — 1 we have A a £ K . So reduction and modular inversion 
in |Lcmma l| are just evalutation of polynomials and inversion in K, respectively. 
Hence 



A a = (f(w) J] (1 - aw) 

= (f(a- 1 ) n a-o)' 1 1 



nfia- 1 ) 



As A(0) — J2 a ™= a A a the lemma follows. □ 
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Lemma 2 is provided here mostly as a starting point for further investiga- 



tions. In the version of the algorithm explained below we will decompose 1 — aw" 
into linear factors and treat all the poles individually. 

Another important consequence of |Lcmma 1 1 is that the type of partial frac- 
tion decompositions Q used to compute vector partition functions can be com- 
puted "for formal 6" in the following sense: 

Proposition 1. Let K be a field and w a variable. Fix a G K x , m S N, and 
f(w) £ K[w] such that f(0) ^ and /(a -1 ) + 1 0. For each b £ N consider the 
partial fraction decomposition 

1 Ax(b;w) A 2 (b;w) B(b;w) 



(1 - aw) m f(w)w b (1 - aw) m f(w) w b 
Then Ai(b;w) is given by an effectively computable expression in b. 



Proof. By Lemma 1 Ai(b; w) is the inverse of f(w)w b modulo (1 — aw) m . Let 
w be the local coordinate at a , i.e., w = a -1 + u>, and carry out all further 
computations in K[w]/(w m ). Then the inversion problem is reduced to invert 
/(a" 1 + w)(a- 1 + w) b = /(a" 1 + w) a - b (l + aw) b . Let g(w),h(w) £ K[w]/{w m ) 
be the inverses of /(a -1 + w) and 1 + aw, respectively. Then g(w) can be 
computed explicitly. The inverse of /(a -1 + w)(a~ 1 + w) b is g(w)a b h(w) b . Note 
the general formula 

JL -^fi + b - l \^-i , , & (&+i) 

By ^ we have 

Hence we obtain explicit expressions for the coefficients of h(w) b , hence for 
g(w)a b h(w) b . Substituting w — w — a" 1 we obtain Ai(b; w). □ 

Example 1. Suppose we want to compute the first numerator of the partial 
fraction decomposition 

1 _ A(b;w) B(b;w) 



(l-w) 2 w b (l~w) 2 



By Lemma 1 we have to invert w b modulo (1 — w) 2 . We do this as indicated 
in the proof of proposition 1| Restated in the local coordinate at 1, we have to 
invert (1 + w) b modulo w 2 . This inverse is (1 — w) b = 1 — bw, a special case of 
([7| . Substituting w = w — 1 we obtain the result 

A(b; w) = l- b(w - 1) = b + 1 - bw. 

Note that the constant term is A(b; 0) = 6+1. We have thus computed our first 
partition function, namely y>i,i(6) = 6 + 1. 
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Example 2. For a more complex example, consider the partial fraction decom- 
position 



1 _ 1 

(1 - w 2 ){l -w 4 )w b ~ (1 - iw)(l + iw)(l + w) 2 (l - w) 2 w b 

A 1 (b) , A 2 (b) , A 3 (b;w) , A 4 (6;u;) , 5(6; «;) 



1 



1 + iw (1 + w) 2 



By direct application of 



we obtain A\(b) 



To 



and A 2 ( b) = l -g 

We have to 



Lemma 1 

determine A 3 (b;w) we proceed as in the proof of Proposition 1 
invert (1 — iw)(l + iw)(l — w) 2 w b modulo (1 + wj T . In the local coordinate 
w at —1, the first factor is 1 — i(— 1 + w) = (1 + j)(l — yxjw), so its inverse 
is + ixj^)- The second and third factor can be inverted similarly. The 

crucial step is to invert (—1 + w) b = (— l) b (l — w) b . This results in (— l) _f, (l + 
w) = (— l)~ b (l + bw). One would proceed by multiplying the four inverses 
modulo w 2 and changing the coordinate back to w, thus obtaining A^{b\w). 
The remaining numerator ^4 4 (6;w) can be computed similarly. 

To carry out the complete calculation by hand would be a considerable 
amount of work. Yet I hope to have convinced you by these indications that 



the proof of Proposition 1 does indeed furnish an algorithm to compute formal 
expressions for the numerators. 



4 Constant terms and partial fraction decompo- 
sition 

We continue by making precise in which sense partial fraction decompositions 
can be used to compute the constant term of a rational function. The first step 
is to define the rings and fields involved in the computations. We consider the 
generating function Ja(z) given in (JT|) as a function in the last variable z m . If 
the last row of A only contains the numbers and 1, all poles of Ja{z) lie in 
C(zi, . . . , z m —i). In the general case they are only contained in an algebraic 
extension. This gives rise to the following definition. 

Definition 3. For a field K, let K(('^/w)) denote the field of formal series 

oo 

f(w) = J2 C k/NW k ' N 
k—M 

for some M € Z, N € N with coefficients c^/jy € K. Similarly let if[ c v/w| 

denote the ring of formal finite sums X)fc!=M c k/NW k ^ N for some M,M' e Z, 
N > 0. Let K(^/w) denote the quotient field of K[rfw\. 

Beware that K((^/w)) is not complete in the natural topology. For example, 
the series J2T=2 wk+l; does not converge. 
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Given a rational function f(w) £ C(w) in a complex variable w, there is 
a Laurent series expansion for each open annulus not containing any pole of 
f(w). We define the "standard" expansion to be the one on the annulus {z £ 
C : < \z\ < e} for some e > 0. If we deal with rational functions over arbitrary 
fields, or even functions involving rational exponents as defined above, there is 
no notion of convergence. Yet we can define the standard expansion as follows: 

Definition 4. Let K be a field and w a formal variable. The standard expansion 
is the embedding K(°^/w) — > K(fflw)) over K which maps w q to w q for all 
q € Q. The constant coefficient of f(w) £ K(^/w) is the coefficient of w° in the 
standard expansion. It denoted by const™ f(w) £ K. 

Using this definition we can state the lemma underlying the partial frac- 
tion method precisely. With regard to the algorithmic application, we use the 
maximal abelian extension Q ab of Q as a ground field instead of C. 

Lemma 3. Consider a vector partition junction ipA- Suppose that A £ Nq™'' 4 '' 
has nonnegative coefficients. Let /a(^) £ Q ab ( c ^zT, . . . : 'y/zm) be the generating 
function of ifA as i n 

(JlJ. Then 

tfA{b) = const Zl • • • const Zm f A{z)z~ h . 

for all b £ Q m . 

We will write const z = const Zl • • • const Zm for this type of iterated constant 
coefficients. The basic principle to successfully apply the partial fraction method 
is to alternate working in fields of type Ki^/w) and fields of rational functions 
K{w). Namely, partial fraction decomposition is applied to rational functions, 
but produces expressions with rational exponents. The following trivial lemma 
allows to clear out the rational exponents for the next iteration. 

Lemma 4. Let f(w) — J2q£Q c g w<! ^ -^(Cv^)) an d r e Q- Let f(w r ) = 
J2q£Q c q w rq . Then const™ f(w) = const™ f(w r ). 

We are now in position to make the statement underlying the partial fraction 
method precise. For q £ Q, let e(q) — e 27rtq . To refer to sets of variables, we 
use the notations z — (z%, . . . , z m ), z' = (z\, . . . , z m —i) and w = z m . Similarly, 
c' k denotes the vector obtained from Ck by omitting the last component, so 

Proposition 2. Let A £ N ( ™' d} , b £ Z m , and q £ Q d . Suppose that the 
columns c\,...,Cd of A are pairwise linearly independent, that a m k > for 
k = 1, . . . , d, and that b m > 0. Consider f(z) = fA,b.q( z ) £ Q ab (z) of the form 

Then one can effectively compute c n £ Q ab , A^ £ N^" -1 *^, 6 (n) £ Z™" 1 , and 

q (n) e Qd j or n = ^ _ _ _ ) N such that 

const z f(z) = c\ const z ' fi(z') + ■ ■ ■ + cn const z ' /iv(z ) 
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for f n (z') = /AC«),&(»), g (»)(z')- 

Note that among the restrictions only the one about the linear independence 
of the columns is severe: If b m < we know immediately that const z fA,b,q(z) = 
0. The factor corresponding to a possible k with a m k = is contained in Q ab (V), 
hence just a scalar. 

Proof. Consider f(z) = f(z';w) as a function of w. All poles of f(z';w) are 
contained in Q ab (^z 7 ). More precisely, the roots of 1 — e(qk)z Ck are w = 

e ( a"^ 1 ) z ' Ck ^ amk for I = 0, . . . , a m k — 1- Hence this factor decomposes as 

This gives rise to the partial fraction decomposition 

, S(6,z';u;) ^ a ^ 1 A fe ;(6;z') 

/a,&,«(*;«0 = — is — + 2^ 2^ 



Here deg^, B(b, z'; w) < b m , so the corresponding summand does not contribute 
to const w fA,b,q(z'] w). As the remaining poles are simple, the numerators 
A k j(b; z') do not depend on w. Hence 



d Q m i.-1 



const w fA,b,g( z> ' w ) = ^ ^2 A k:l (b;z'). 

k=0 1=0 

It follows that 

d a TO n,- l 

const z fA,b, q {z'\ w) = const 2 / const™ fA,b,q{z'\ w ) = X! X! const ^' ^)m(&; z')- 
Our strategy to construct fi(z'), . . . , /at(V) is hence to transform the summands 



Ak.i(b; z') suitably. By Lemma 1 they are given by 



A kjl (b;z') = l / l[(l-e(q~ k )z a «) 

I k^k 



11 1 

l^l \ 



After the substitution, the first factor in the denominator becomes 

n( ( . a mk(i- qk )\ ,'- 
1 - e I q~ k + ] z 



k^k 



M ink 
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the second factor becomes 

' b m (l - q k ) 

&mk 

and the last factor becomes 




We want to construct c\, . . . , e/v and f±(z'), . . . , f^(z') from Ak.i(b, z'). To sim- 
plify notation, change the indexing to n — (k,l). The idea is to choose 

1 

Ck,l 



&rnk 



and 



A,lCO = —A k .i(b- z') 



1 - e ?S 



a mk{l-q k ) 
^mk 



This is not quite possible, because b^ k ^ = (b m c' k — a m kb') / a m k contains rational 
numbers as opposed to integers, and the columns (a m j,c~ — a rn k c 'k)/ a mk of A^ k,} ^ 
contain rational numbers as opposed to nonnegative integers. By a multivariate 
version of 



Lemma 4 



we can substitute Zj by z^ mk without changing the con- 
stant coefficient, thereby obtaining integers in both places. To to make A^ k ^ 



nonnegative, use the formula 



7~ b 1 - Z b 



(8) 



as necessary. The sign will amount to a change of Ck,u the numerator to a 
change of U k ' l \ □ 



The most important feature of Proposition 2 is that the algorithm does 
not branch on a condition depending on b. Hence, although stated for specific 
b G Z", it can actually compute formal expressions in b. 

Example 3. Let A = (3 })• Then <pA(a,b) = const /a{z, w)z~ a w~ b , and 
f A (z,w)z~ a w~ h 



(1 — zw 3 )(l — zw)z a w b 



This is a trivial example, included here only for the sake of illustration. In the 
first step of the iteration we compute the constant term with respect to w. It 
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b < then (Pa(o>, b) — 0, so suppose b > 0. As z a is a scalar now, we start with 
the following partial fraction decomposition specializing ^ : 

1 = Arjb-w) | A 2 (b;w) | B(b;w) 

(1 — zw 3 )(l — zw)w b 1 — zw 3 1 — zw w b 

In order to apply [Lemma 2 we have to ensure that z is a third power. So let 
z 1 / 3 be a third root of z, and use if = C(z 1 / 3 ) as the field of scalars. Then 

const™ - ~-r —r = Ax{b; 0) + A 2 (b; 0) 

(1 — zw 6 ){l — zw)w° 

— S(i-zw)w b ,zvV + 5 f (i- 2U ,3) UI !)^(l). 
Let ( be a primitive third root of unity. Then 

S(l- zw)wb JS) = - ^ (1 _ z 2/ 3)z - b / 3 + (1 _ Cz 2/3) C 6 z -b/3 

1 \ 



(1 -(2 z 2/3^2b z -b/3 J ' 
B(l- zv fi) w *,zW (l_ z -2) z -b (l-z 2 ) Z - b - 2 ' 

Hence ip^(a,b) is the constant coefficient of the univariate function 
if I 1 



3 V (1 — z2/3) z a-6/3 (! _ ^ 2 2/3)^6 z a-b/3 

1 



-6-2 ' 



(1 - C2^2/3^26 za -6/3 J (J _ z 2) z a 

In order to eliminate the fractional powers, we substitute z by z 3 according to 
ILcmma"4l in the first three summands and obtain 

, n if 1 1 

<Pa{0; b) = const 



3 1(1 - z 2 )z 3a - b (l-(z 2 )( b z 3a - b 



(l-( 2 z 2 )( 2b z 3a - b J {l-z 2 )z 

From this we would compute a quasi-polynomial expression for ^(a, b) as a 
second step. 

Example 4. Let A = (J ? i )• Then ip^ is Kostant's partition function associ- 
ated with a root system of type A 2 . A prototype implementation of the partial 
fraction method produces, given A, the following output: 



a + 1 if a > If a - 6 i£ a - 6 > 1 
^(a, 6) = ^ 1 otherwise J 1 otherwise 



if 6> 
otherwise. 
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Note that the traces of the algorithm are still visible in the output: The output 
consists of nested case statements, where the second alternative is always 0. This 
corresponds to to the fact that if b m < then const z f(z)z~ b = 0. The outer 
case statement corresponds to the first step, the partial fraction decomposition 
with respect to z 2 . Hence the condition is simply b > 0. As there are two 
nonzero components in the second row of A, the first iteration yiels a sum of 
two rational functions in z\. These correspond to the two terms of the difference. 
Note that for the second term there is a change of sign, and the condition is that 
a certain function of a and b must be positive. This is because of an application 
of ^ in the first step, as explained in the end of the proof of Proposition 2 
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